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ABSTRACT 



We examine the theoretical relationship between fio and substructure in galaxy 
clusters which are formed by the collapse of high density peaks in a gaussian random 
field. The radial mass distributions of the clusters are computed from the spherical 
accretion model using the adiabatic approximation following Ryden & Gunn. For a 
cluster of mass, M{r,t), we compute the quantity AM/ M at a cosmic time t and 
within a radius r, where AM is the accreted mass and M is the average mass of 
the cluster during the previous relaxation time, which is computed individually for 
each cluster. For a real cluster in three dimensions we argue that AM/ M should be 
strongly correlated with the low order multipole ratios, (I)^*/'^™*' ^^le potential 
due to matter interior to r. Because our analysis is restricted to considering only the 
low order moments in the gravitational potential, the uncertainty associated with the 
survival time of substructure is substantially reduced in relation to previous theoretical 
studies of the "frequency of substructure" in clusters. 

We study the dependence of AAI/ M on radius, mass, fioi Aq = 1 — r^O: redshift, 
and relaxation timescale in universes with Cold Dark Matter (CDM) and power-law 
power spectra. The strongest dependence on f2o (Aq = 0) occurs at z = where 

AM/ M cx for relaxation times ~ 1 — 2 crossing times and only very weakly 

depends on mass and radius. The fractional accreted mass in CDM models with 
^^0 + Aq = 1 depends very weakly on fig and has a magnitude similar to the 
ilo = 1 value. AM/ M evolves more rapidly with redshift in low-density universes and 
decreases significantly with radius for fig = 1 models for z ^ 0.5. We discuss how to 
optimize constraints on fio and Aq using cluster morphologies. 

It is shown that the expected correlation between AM/ M and $»"*/$«"* extends 
to the two-dimensional multipole ratios, ^'™*/^o"*' which are well defined observables 
of the cluster density distribution. We describe how N-body simulations can quantify 
this correlation and thus allow AM/ M to be measured directly from observations of 
cluster morphologies. 

Key words: galaxies: clusters: general - galaxies: evolution - galaxies: structure - 
cosmology: theory - X-rays: galaxies. 



1 INTRODUCTION 

Clusters of galaxies have proven to be useful laboratories for 
cosmological studies and have in particular yielded interest- 
ing measurements of f2o, the present value of the average 
mass (energy) density of the universe expressed in terms 
of the critical density required for closure. The methods to 
measure fio with clusters each have their own advantages 
and disadvantages (see Dekel, Burstein, & White 1997). Per- 
haps the most familiar method is to assume the mass-to- light 
ratio in clusters is representative of the universe as a whole 
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( ICarlb erg, Yee, fc Ellingson 1997|) . This method is popular 
because of its conceptual simplicity and relative ease to im- 
plement, but the basic assumption of this method may be 
unjustified. Another promising method uses the abundances 
of clusters (e.g. White, Efstathiou, & Frenk 1993) to derive 
a relationship between f2o and erg, the rms fluctuations in 
mass within spheres of radius ?>h~^ Mpc. This method has 
the advantage that the theoretical comparison can be use- 
fully achieved by relatively simple and computationally in- 
expensive semi-analytical methods. One disadvantage of this 
method is that cluster masses are usually determined from 
X-ray temperatures, the relationship of which is not entirely 
certain. Another disadvantage is the uncertainty in due 
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the un certa inty in the bias factor for clusters (Mo, Jing, & 
White 



997) 



The detailed structure of clusters provides a comple- 
mentary constraint on Q.q. Although the radial density pro- 
files of clusters appear to be very insensitive to f2o when 
scaled in terms of their virial radii (Navarro, Frenk, & White 
1997), the presence of substructure in the mass distribu- 
tion appears to be a quite sensitive diagnostic of f2o as sug- 
gested in the pioneering study by Richstone, Loeb, & Turner 
(1992). Advantages of this method include the direct influ- 
ence of fio on cluster morphologies (i.e. growth of structure 
regulated by f2o), and cluster morphologies are straightfor- 
ward to quantify and compute observationally (e.g. Buote 
& Tsai 1995). The principal disadvantage is that a proper 
theoretical model requires computationally expensive high- 
resolution N-body simulations of a large number of clusters 
sufficient to adequately sample the variety of spatial mor- 
phologies of a cluster population (Jing et al. 1995; Buote & 
Xu 1997; Thomas et al. 1997). 

To date the only study that has compared the requi- 
site large N-body simulations to observations is Buote & 
Xu (1997; hereafter BX). Using dissipationless simulations 
of scale-free and Cold Dark Matter (CDM) models (e.g. Os- 
triker 1993) and X-ray data of clusters from ROSAT, BX de- 
termined that the spatial morphologies of galaxy clusters (1) 
decouple the influence of f2o and the power spectrum when 
quantified in terms of statistics based on ratios of multipoles 
of the projected gravitational potential (Buote & Tsai 1995) , 
and (2) favor a low matter density (fio ^ 0.3). However, a 
rigorous determination of fio from X-ray cluster morpholo- 
gies requires dissipational simulations of a statistically large 
number of clusters, at present a daunting computational ex- 
penditure even for only one set of model parameters^. To 
facilitate this process, less computationally intensive semi- 
analytic methods could be used to rapidly identify the op- 
timum parameters (e.g. redshift, mass range, radius) of a 
cluster sample which maximize the ability to determine f2o 
and possibly a cosmological constant via cluster morpholo- 
gies. 

A re-examination of semi-analytica^ models of the rela- 
tionship between cluster morphologies and fio is important 
in its own right. As it is to be expected given the nature of 
the problem, semi-analytic approaches to cluster morpholo- 
gies have been highly idealized. The previous studies do not 
actually construct model cluster mass distributions having 
substructure, but either adopt a spherical collapse model 



for co mputing a d istribution of cluster collapse times (Rich- 
stone e: al. 1992) or use simplified models of the merging 



histories of clusters to compute the distributions of times 
since a substantial merger event (Kauffmann & White 1993; 
Lacey & Cole 1993). Even if we accept such simplifications 
as necessary evils to be traded for intuitive guidance and 
increased computational speed, these models are of limited 
usefulness because they only compute the ambiguous "fre- 
quency of substructure" in clusters which is very sensitive 



t Dissipationless simulations should be adequate for mass maps 
of clusters derived from gravitational lensing when such maps for 
observed clusters become widely available. 

i For our discussion we broadly define semi-analytical to be tech- 
niques other than three-dimensional N-body simulations. 



to the uncertain survival time of substructure (Kauffmann 
& White 1993; Lacey & Cole 1993). 

In this paper we present an intuitive model for substruc- 
ture, based on the spherical accretion model, that is closely 
related to quantitative observables of cluster morphology. 
We describe how N-body simulations can quantify the re- 
lationship between our model and the observables and thus 
allow a direct measurement of fio- However, for our present 
study, we use this model to (1) provide physical insight into 
the evolution of cluster morphology with redshift as a func- 
tion of f2o and Ao = 1 — fio and to (2) assist the designing 
of future observing programs for the detailed comparison of 
these observational samples to clusters produced in N-body 
simulations. We then outline a procedure to optimize de- 
termination of S7o and Ao using the observables of cluster 
morphology that are closely related to our model. 

The paper is organized as follows. In §^ we motivate our 
model and show its relationship to quantitative observables 
of cluster morphology. The details of the spherical accretion 
model and the cosmological framework are discussed in §^ 
and §^ respectively. In §^ we analyze the model for different 
values of f2o, Ao, and z for CDM and power-law power spec- 
tra. Finally, we discuss these results and identify the best 
ways to constrain flo and Ao observationally in «J. 



2 MOTIVATION 
2.1 Theory 

In a standard Friedmann-Robertson- Walker universe with 
S7o < 1 and Ao = 0, the linear growth of density fluctuations 
becomes strongly suppressed when the curvature term in the 
Friedmann equation exceeds the matter term. The redshift 
delineating this transition from an Einstein - de Sitter phase 
to one of free expansion is then 1+ztrans = — i.e. when 
the matter density fl{zt,-ans) ~ 0.5. (e.g. §11B in Peebles 
1980). Hence, if fio ^ 1, then objects formed a long time in 
the past relative to universes where f2o ~ 1, and thus clusters 
in low-density universes should be, on average, more relaxed 
than clusters in universes with Slo ~ 1- To apply this idea 
to real clusters it is necessary to specify what it means for 
clusters to be "more relaxed". 

The most general stable self-gravitating non-rotating 
equilibrium conflguration is the triaxial ellipsoid. Features 
like substructure that break this symmetry provide a mea- 
sure of the departure of a cluster from a virialized state. The 
relaxation rate determines how rapidly substructure and 
other non-ellipsoidal features are erased in order to bring 
the system to equilibrium. Hence, the amount of substruc- 
ture in a cluster at a particular epoch, or rather the degree 
of non-ellipsoidal symmetry, is approximately determined by 
the amount of mass accreted over the timescale associated 
with the relaxation rate. 

This relaxation timescale depends on the mass distri- 
butions of the cluster and the subclumps (e.g. White & 
Rees 1978; Binney & Tremaine 1987), and is typically of 
order the crossing time of the system for substructures and 
non-ellipsoidal distortions comprising a substantial fraction 
(^ 20%) of the cluster massPI. Since our interest lies in how 



Mergers of structures of these sizes with clusters are typical for 
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much a cluster departs from a virialized state, we must con- 
sider how much a cluster evolves morphologically over the 
relaxation timescale, not just the amount of substructure 
present at a particular time. 

If gravity drives the dynamical evolution of clusters, 
then we would expect the morphological evolution to be 
similar for clusters of different masses (approximately self- 
similar); i.e. the degree of morphological evolution should 
be proportional to the amount of accreted mass and in- 
versely proportional to the total cluster mass. We write this 
fractional amount of mass accreted over a crossing time as 
AM/ M, where AM is the accreted mass and M is the av- 
erage mass of the cluster within a radius r. We may write an 
equivalent expression in terms of the multipole moments of 
the gravitational potential due to matter interior to a radius 



AM 

M 



A$y' 

—int ' 
$0 



(1) 



where $o"' — —GM{< r) /r is the monopole term. Of course, 
substructure and other non-ellipsoidal distortions have im- 
portant contributions to other moments of the potential, and 
we expect that the increase in the monopole will be strongly 
correlated with increases in the next few multipoles. (Indeed, 
this is the basis of our hypothesis that substructure is re- 
lated to the amount of accreted mass over the relaxation 
timescale.) Hence, let us consider the increase in <|>'"'j ex- 
panded in terms of all of its multipoles (and suppressing the 
azimuthal terms). 



int 
I 



(2) 



where "!>;"' is the Ith term in the multipole expansion at 
time t of interest, and q; characterizes the increase in 
over the previous crossing time; i.e. a; = if the particular 
multipole does not increase, and a; = 1 if the multipole term 
is zero at time t — tcross and non-zero at t. For our discussion 
we will consider a spherical region positioned at the center of 
massfj] and thus we do not require the azimuthal terms in the 
multipole expansion - although our general argument does 
not depend sensitively on the shape of the aperture used 
to compute the moments. Since the I > terms are not 
spherically symmetric, we average over the spherical 

surface of radius r and denote it by {$"''). Actually we 
average ($j"')2 because = for « > 0. 

Because we cannot observe a cluster for the duration of 
a crossing time, let us approximate equation at the time 
t of interest. First, in the denominator we set ~ 
which is accurate to within a factor of 2, and set « 
^o"', since the monopole term dominates the higher order 
terms for all reasonable cases. For a relaxed, smooth cluster 
the / > ratios {($r')^)/((*o"')^) are substantially smaller 
than for a cluster with subclumps that are a considerable 

(^ 20%) fraction of the total mass (Buote & Tsai 1995)[[Jj 



clusters formed in N-body simulations (e.g. Tormen, Bouchet, & 
White 1997). 

^ The center of mass can be defined within the sphere of radius 
r through iteration. 

I Buote & Tsai actually discuss circularly averaged moments of 



As a result, we may approximate equation (g) by setting 
ai — 1 for the I > terms Making these substitutions 
we arrive at, 



/AM\ 



+ 



((<I.-')2) ■ 

l>0 >^ u ' ' 



(3) 



which states that the fractional increase in the rms spher- 
ically averaged potential over a crossing time is approxi- 
mately the fractional increase in the mass and the increases 
in the ratios of the rms spherically averaged higher order 
multipoles to the monopole added in quadrature. Violent 
relaxation ( Lynden-Bell 1967 ) is the key process driving the 
elimination of large potential fluctuations. It operates on a 
timescale of ~ 1 — 2 crossing times and proceeds indepen- 
dently of the masses of the constituents in accordance with 
our picture of the morphological evolution of clusters stated 
above. 

Hence, equation (^) is a definition of the dynamical 
state of a cluster and could have been the starting point 
of our discussion^ By the nature of the multipole expan- 
sion, we see that only the low order moments are important 
for characterizing the dynamical state, with the monopole 
term being most important^. As a result, the possible long- 
term survival of the dense, compact cores of subclusters, 
which contribute mostly to the higher order moments, does 
not affect the shorter relaxation timescale of the low order 
moments most relevant to the dynamical state. (The impli- 
cations of this property are discussed in §§0 and 0.) 



2.2 Observational Consequences 

The fractional change in the monopole, AM/ M, cannot 
be directly observed over a relaxation time. However, our 
premise that the amount of accreted mass over the pre- 
vious relaxation timescale determines the amount of sub- 
structure (or non-ellipsoidal distortions) requires that the 
monopole change be strongly correlated with the change in 
the low order terms, which are approximately the ratios, 
{(*r')^)/{(*o"')^), defined at the epoch under considera- 
tion. What is the nature of this correlation and, more im- 
portantly, how may AM / M he inferred from cluster obser- 
vations? 



the projected potential which replace the cluster mass density 
with toy models of the X-ray surface brightness of clusters. Nev- 
ertheless, the qualitative features apply when the mass is used as 
shown §4 of Tsai & Buote ( |l99e| ). 

** This approximation implies that ((A<I>™')2)/(('I>'"')2) for 
virialized clusters because the even-l moments, though small in 
relation to clusters with considerable substructure, are non-zero 
in general for ellipsoidal masses. 

tt For alternative indicators of the d ynamical state of a cluster 
see Zaroubi, Nairn. &; Hoffman ( |l996| ) and Natarajan, Hjorth, Sz 
van Kampen (|l997|). 



M Our formulation is a twist on the argument of Richstone et al. 
(1992) who computed the fraction of present-day clusters which 
"collapsed" over the previous relaxation time. This fraction is 
then interpreted as the "frequency of substructure" in the present 
cluster population. In our case, for an individual cluster the frac- 
tional accreted mass over the relaxation time approximately de- 
termines the importance of substructure in the cluster and thus 
its dynamical state. 
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For observations of clusters it is convenient to work 
with the projected potential, and its circularly aver- 

aged (over radius R) multipole ratios, {(vE'm*)^)/{(*o"'')^)- 
For the m > terms it is unclear how to relate a mea- 
sured ((*™')^)/{(*o"')^) to the multipole ratios in three- 
dimensions. However, the fractional change of the monopoles 
in 2-D and 3-D are comparable if the fractional accreted 
masses within the 2-D and 3-D regions are comparable. That 
is, let us again consider our prototype case where the 3-D 
region is a sphere of radius r defined at its center of mass. 
The corresponding 2-D region is the circle of radius R = r 
resulting from the projection along its symmetry axis of the 
cylinder which encloses the sphere at its center. So long as 
the extra accreted mass in the portion of the cylinder outside 
the sphere does not contribute much to the total fractional 
accreted mass of the cylinder, then the 2-D and 3-D frac- 
tional monopole changes will be similar. 

Hence, A^j,"'/ ^ for projected radii 

containing a sizeable fraction of the total cluster mass. (An 
appropriate minimum radius can be determined from N- 
body simulations.) If we now consider the 2-D analogue of 
equation (^), then the correspondence of the 2-D and 3- 
D monopole changes coupled with the arguments made in 
the preceding section indicate that (AM/ M)^ (3-D) should 
be strongly correlated with {(*m')^)/((4'o"*)^)i which are 
themselves easily computable observables of the cluster den- 
sity distribution (Buote & Tsai 1995). For observations of 
an individual cluster this correspondence assumes that pro- 
jection does not smear out the structure; e.g. axisymmetric 
clusters are viewed edge-on. Later in this section we discuss 
how to deal with projection effects^. 

Three dimensional N-body simulations can explic- 
itly quantify the relationship between AM/ M and 
((*m*)^)/{(*o"')^)- There is reason to believe that this re- 
lationship will be very tight for small m and mostly inde- 
pendent of fio. For clusters observed with ROSAT (Buote 
& Tsai 1996) and those formed in both dissipationless 
(BX; Thomas et al. 1997) and dissipational N-body simula- 
tions ( [Tsai fc Buote 199e| ), the ratios, ((*j;^*)^)/{(*j,"')^) , 
are strongly correlated for the first few m - in fact the 
logarithms of these ratios are essentially proportional to 
each other P|. We henceforth adopt the simplified notation 
P™/Po = {(*^*)^) Ailp"*)^) used by these other studies. 

Tsai & Buote (1996) studied these correlations in detail 



for both the X-ray surface brightness maps and the projected 
masses of six clusters formed in the N-body / hydro dvnam- 
ical simulation of Navarro, Frenk, & White (1995). These 
correlations can be understood as tracks followed by clus- 
ters as they evolve via mergers and subsequent relaxation 
from infancy to quasi-virialized states (see Figure 6 of Tsai 
& Buote and §5.1 of Buote & Tsai 1996). For the larger val- 
ues of Pm/Po (corresponding to the least relaxed clusters) 



the m = 2, 3, 4 ratios are tightly correlated and are very 
nearly proportional to each other. 

The scatter in these correlations is largest for the 
smallest values of Pm/Po as a result of projection effects. 
The smallest values of Pm/Po correspond to either relaxed, 
single-component clusters or to clusters with substructure 
smeared out by the act of projecting along the line of sight; 
e.g. a bimodal cluster being viewed along the merger axis. 
Thus, to minimize the projection effects the correlations be- 
tween the low order Pm/Po (particularly for logj^Q P2/P0 ver- 
sus logj^Q P4/P0) can be defined accurately for the larger val- 
ues of Pm/Po where the scatter is low (see, e.g.. Figure 3 of 
Buote & Tsai 1996, Figure 4 of BX, and Figure 9 of Thomas 
et al. 1997.) Moreover, the approximation a; = 1 in equation 
(^) is most rigorously satisfied for these largest multipole ra- 
tios. 

The shapes of these correlations do not show siginfi- 
cant differences as a func tion o f fio or Ao = 1 — f2o (BX), 
although Thomas et al. (1997) note that the logjQP2/Po 
versus logj^g P-i/Po correlation is tighter for the Ao = mod- 
els. Since these correlations are expected to apply as well to 
AM/M (3-D), we may write. 



(4) 



where Cm and dm are the linear coefficients which should be 
largely independent of f2o and, for best accuracy, should be 
determined from fits to the largest values of log^g Pm/Po- By 
quantifying the Cm and dm with N-body simulations we may 
infer AM/ M from observations of Pm/Po for a large sample 

of clusters and thus measure the value of Qof^ (Another ad- 
vantage of writing the expected correlations in terms of log- 
arithms is that BX found that the means of the log^g Pm/Po 
(m — 2, 3, 4) distributions of N-body clusters are infiuenced 
primarily by Slo whereas the variances are affected mostly 
by the power spectrum of density ffuctuations; i.e. the inffu- 
ence of the power spectrum may be largely circumvented by 
examining the means of the logarithmic distributions.) 

Measurements of Pm/Po for clusters may be obtained 
from either gravitational lensing or X-ra y data (for a discus- 
sion see §2 of Buote & Tsai 1995)* * *. Since gravitational 
lensing measures the projected mass directly, the coefficients 
Cm and dm can be quantified by dissipationless simulations 
which is an advantage over X-ray data. However, the pro- 
jection of the mass in simulated clusters appears to lead to 
smaller measured variations in Pm/Po than does the X-ray 
surface brightness (see §4 of Tsai & Buote 1996), and thus 
the larger dynamic range afforded by X-ray maps probably 
indicates a greater sensitivity to flo- This larger dynamic 
range occurs because for unrelaxed clusters the X-ray emit- 
ting gas approximately traces the mass but its emission is 
more responsive to density ffuctuations since it goes as p'^as- 
In the other limit, when the gas is in hydrostatic equilibrium 
the emission traces the potential for arbitrary temperature 



§§ We show in § ^.2.l| that the {{■^'t^'^f )/{{%"'' f) obtained from 
the N-body simulations of BX agree reasonably well with our 
calculations for AM/ M. 

Again, although most of the references listed have focused on 
moments of the X-ray surface brightness, Tsai & Buote (1996) 
showed that taking moments of the mass gives very consistent 
results. 



Ill A minimum of about 20 clusters appears to be required to 
distinguish between CDM models with Qq = 1 and Co = 0.3 
(Thomas et al. 1997). 

" * Another route to observations is via cooling flows (e.g. 
Fabian 1994), where t he cooling-flow rat e should be anti- 
correlated with AM/ M ( Buote fc Tsai 199{ ) 
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gradient (§3.1 of Buote & Canizares 1994; §5.1 of Buote & 
Canizares 1996) and is thus smoother and rounder than the 
mass (see §2.3.1 of BX for a detailed discussion). 



2.3 Present Application 

Presenting a constraint on Qq by obtaining the coefficients 
Cm and dm (equation ^ from N-body simulations for com- 
parison to observations is beyond the scope of the present 
paper. Instead we prefer to focus on the differences be- 
tween CDM models with different Qo since the differences 
in AM/ M can be related to differences in Pm/Po with- 
out knowledge of the as yet unknown coefficients dm and 
by taking c™ « 1 as suggested by the logj^Q Pm/Po correla- 
tions. For this task we can compute the quantity AM/ M 
using a spherical model for a cluster which has the signifi- 
cant advantages of being simpler conceptually and consider- 
ably less intensive computationally than three-dimensional 
N-body simulations. In this way we can gain physical insight 
into the morphological evolution of clusters in different cos- 
mologies and rapidly explore the interesting parameter space 
to aid analysis of future three dimensional simulations and 
observations. 

Of course, a spherical model is only an idealization, but 
the radial density profiles seen in three dimensional N-body 
simulations agree well with simple spherical accretion mod- 
els (e.g. Navarro, Frenk, & White 1995; Anninos & Norman 
1996; also see §Bj| of this paper). By using a spherical model 
we in essence compute a mean value of AM/ M, whereas in 
a more realistic model the non-spherical nature of the merg- 
ing process will induce fiuctuations in AM/ M depending 
on how the mass is distributed in clumps during a merger. 
However, for studying the dependence of AM/ M on f^o, 
the mean value should be sufficient since BX showed that 
it is primarily Qo which determines the means of the distri- 
butions of Pm/Po for clusters formed in three-dimensional 
N-body simulations. We therefore adopt a spherical model 
for our present investigation. 



3 COMPUTATION OF AM/M 
3.1 Peaks formalism 

It is our intention to work within the standard framework 
wherein small density inhomogeneities in the early universe 
grow by gravitational instability into the galaxy clusters that 
we see today. We restrict our discussion to the case where 
these initial density fluctuations are described by a gaus- 
sian random fleld, and, in particular, that the sites where 
clusters form are determined by high peaks in this pertur- 
bation field (Bardeen et al. 1986, hereafter BBKS). With 
these assumptions the average initial density profile around 
such high peaks is specified by the power spectrum of density 
fluctuations, P{k). (For similar presentations of the follow- 
ing see Ryden & Gunn 1987, Ryden 1988, Hoffman 1988, 
and especially Lilje & Lahav 1991.) 



The spherical accretion model (see §3.2) only depends 



on the density through the initial cumulative density con- 
trast, Si{x), defined at comoving position x at initial time ti. 
The initial cumulative density contrast around a high peak 
in a gaussian random field is (BBKS), 



1 — 7 / cro 1—7^ (72 



(5) 



where 7 = cri/((T2cro) and the spectral moments are defined 
by, 



<T- = ^ I P(k)k'^-'^dk, 



(6) 



with the special case (jg = ^(0); the cumulative correlation 
function and its cumulative Laplacian are. 



^ " 2^7 

VH, = ^ J P{k)k'^w{kx)dk, 

(note: V^^ ~ where the top-hat filter is given by, 

3 



'w{x) 



{sinx — X cos x) . 



(7) 
(8) 

(9) 



The dimensionless parameter f specifies the peak height, 
Si{0) = i^cTQ, and 6(7,71^) is given by equation (6.14) of 
BBKS with the asymptotic behavior 6 as u ^ 00. 
Thus, for high peaks 5 ~ v^/ao- The peak only dominates 
the collapse in its vicinity which we take to be defined as 
the region where &i < ^ (AcSi)^, where the variance in &i is 
(BBKS), 

-2 



{A&if = all - 



■r 



n cr2 



27 



0"2 



where. 



1 

2^ 



P{k)w^{kx)k^dk, 



(10) 



(11) 



is the mass variance in spheres of comoving radius x. 

We do not consider secondary perturbations arising 
from non-sphericity and random velocities (BBK S) dis- 
cussed by Ryden & Gunn ( [l987| ) and Ryden ( [l988| ) in the 
context of galaxy-sized halos. Such perturbations cause fiuc- 
tuations in the collapse time of matter around the peak, 
though we do not expect angular momentum to play a dom- 
inant role in the formation of clusters since they are not ob- 
served to rotate appreciably. At any rate, uncertainties of 
this variety that affect the collapse time are related to our 
ignorance of the precise relaxation timescale over which low- 
order potential fluctuations are erased, the uncertainties of 
which we expect to dwarf those of the secondary perturba- 
tions. (We discuss this issue further in the next section.) 

For our calculations we define the initial epoch to be 
at recombination, Zi = 1300, though our results differ neg- 
ligibly for any Zi 50. In practice we compute Si by first 
evaluating So, the cumulative density contrast extrapolated 
to the present using linear theory, by normalizing the lin- 
ear power spectrum at 2; = (see §^). Then Si is computed 
from So assuming linear growth, for which we use th e con - 
venient approximation of Carroll, Press, & Turner (1992) 
for a universe with matter density parameter Q and density 
parameter, A, due to a cosmological constant. 

3.2 Spherical accretion model 

We evolve the initial cluster density distribution determined 
by the peaks formalism to a final state using the spherical 
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accretion model |*~*~*| . It will now prove convenient to work 
in proper coordinates, r = ax, where a(t) denotes the expan- 
sion factor at time t. For arbitrary f2o and a monotonically 
decreasing density profile, the spherical collapse model (e.g. 
Peebles 1980; Padmanabhan 1993; Sahni & Coles 1995) ex- 
actly describes the expansion of spherical shells out to their 
maximum radius ("turn around"). (In the Appendix we give 
the equations for the zero-curvature universe with a cosmo- 
logical constant Aq = 1 — Oq.) After the shell reaches maxi- 
mum expansion, it collapses and crosses the orbits of other 
shells at which point the spherical model no longer accu- 
rately describes the motion of the shell. If one restricts the 
solutions to those that are self-similar, then exact solutions 
may be found (Fillmore & Goldreich 1984; Bertschinger 
1985). For the general case, however, methods must be used 
which specifically treat the interactions between the shells. 
The standard method to do this is with N-body simulations 
(e.g. Sigurdsson, Hernquist, & Quinlan 1995; Thoul & Wein- 
berg 1995). 

Another approach is to consider the growth of a den- 
sity peak by computing the orbit of an infalling shell in 
the potential generated by the previously collapsed matter 
(Gunn 1977; Ryden & Gunn 1987 - also see Blumenthal et 
al. 1986). If the infalling shell induces only a small change in 
this potential, then the shell's orbit may be computed using 
the adiabatic approximation. The accuracy of this approx- 
imation is determined by the number of shells, N , used to 
divide up the mass distribution. Typically, only for the first 
few (innermost) shells does the potential change substan- 
tially. For large A'^ (~ 100) the errors associated with this 
approximation have a negligible effect on shells at radii that 
are of interest to our present study. 

To construct a cluster using the adiabatic method (see 
Ryden & Gunn 1987) one begins with a core mass distribu- 
tion, taken to correspond to the scale over which the power 
spectrum is smoothed (see The details of the shape of 
this core have no tangible effect on the shells at radii of in- 
terest. We consider now the first mass shell orbiting in the 
potential generated by the core mass. The energy per unit 
mass of the shell is an integral of the motion and is equal 
to the potential at the maximum radius of the orbit. More 
generally, we may write for a shell n falling into the poten- 
tial, $,1, generated by the mass of the previous n — 1 shells 
plus the core as. 



P„(r)dr : 



En 



(C) 



dt) 



(12) 
(13) 



where r'^ is the (proper) maximum radius of the nth shell. 
For the first shell we have <E>i = '^c, the potential generated 
by the core mass. The amount of mass a shell contributes 
within a radius interval dr is determined by the fractional 
amount of time that it spends passing through dr during its 
orbit. The amount of time the shell spends within a radius 
interval dr is just dt = dr/vr, and thus the fractional amount 
of time, which is just the probability of finding the shell 
within dr, is, 



* * * For a discussion of some of the caveats associated with this 
type of approach see Lilje & Lahav (1991) and Bernardeau (1994). 



dr/v- 



(14) 



and = \dr/dt\ is the radial speed of shell n, which for the 
first shell is. 



jI = ^2[E^~^,{r) 



(15) 



Hence, the combined mass of the core and the first shell may 
be written generally. 



(16) 



M„(r) = M„_i(r) + Msh / Pn{r)dr 



where M^h is the total mass of the shell and Mo — Mc. 

The next shell is added analogously to the first by com- 
puting its energy, E2, in the potential, $2, generated by 
mass. Ml, and results in a new mass, M2. Now, however, 
the first shell sees a different mass due to the overlapping 
orbit with the second shell. If the potential change is small, 
then the energy of the first shell decreases, but since the 
radial action. 



Vrdr, 



(17) 







is conserved, the effect is to reduce the ma ximum radius of 
the orbit of the first shell (Gunn 1977) * * *. Hence, when the 
second shell is added, jr is recomputed for the first shell in 
the potential generated by the mass M2 . Then r^^ is adjusted 
so that the action equals its initial value. For an adiabatic 
change in the potential, the new value of rm is approxi- 
mately. 



old Jr 



Anew ' 
Jr 



where jr 



Jo 



"dr, with 



(18) 



computed using the 



potential incorporating the new overlapping shell - for the 
case under consideration, the potential generated by the 
mass M2. Equation (^^ can be iterated by computing a 
new jr"™ using r"^™ as the new upper limit in the integral 
until desired precision in rj^*^™ is achieved. Once the radius 
of the first shell has been modified to reflect the addition of 
the second shell, a new value of M2 is computed. 

(At this point one should in principle repeat this proce- 
dure for the first and second shells until r^^ and M2 change 
within desired tolerances. However, corrections of this type 
are negligible when the initial potential change is small, 
which applies as n becomes large. We neglect these higher 
order corrections.) 

This adiabatic addition of shells continues until all of 
the collapsed shells have been incorporated into the aggre- 
gate virialized cluster. In general not all of the mass bound 
to the peak has collapsed at the redshift of interest, say z. 
We define the last collapsed shell to be that which has just 
reached r = at redshift z. We compute this directly by 
integrating dr/vr for each shell rather than estimating the 
collapse time as 2tm- For the shells that are bound to the 

* * * For the first few shells the radii can actually increase be- 
cause the density profile changes drastically in shape with the 
addition of each new shell; i.e. the addition of the first few shells 
is not really adiabatic if the shell mass is comparable to the core 
mass. 
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peak but uncoUapsed at z, we let them fall into the clus- 
ter potential and assign them to their infall radii at z. Note 
that Ryden & Gunn (1987; Ryden 1988) did not treat the 
uncoUapsed matter while Hoffman ( 198^ ) allowed all of the 
matter to collapse. 

After building the cluster mass distribution at the epoch 
of interest, we evaluate the crossing time which we take to 
be, 



,^1 



/2 



(19) 

u 

where R1/2 is the half-mass radius and a is the one- 
dimensional (radial) velocity dispersion computed using the 
virial theorem. Next, we construct the cluster mass distri- 
bution at the earlier time, t — t cross- Finally, we evaluate the 
fractional amount of mass accreted within a radius r over 
the time interval (t — t cross, i), 



AM _ M{r,t)- M{r,t-tcross) 
M ^ l[M{r,t)+A4{r,t-tcross)y 



(20) 



where AM is the mass difference and M is the average 
mass. As we show in the value of AM/ M is quite sen- 
sitive to the time interval {t — tcross,i) used to represent 
the relaxation timescale. We need only consider a relaxation 
timescale that is ~ 1 — 2 crossing times since, as mentioned 
in we are interested only in substructure contributing to 
fluctuations in the low order moments of the gravitational 
potential. 

However, even the relatively modest range of ~ 1 — 2 
crossing times causes sufficiently large changes in AM/ M 
to make further refinements to our spherical accretion model 
unwarranted; e.g. the secondary perturbations due to non- 
sphericity and ang ular momentum discuss ed in §3.2, initial 



peculiar velocities (Bartelmann et al. 1995), or a small (and 
difficult to precisely define) drag force due to substructure 
(Antonuccio-Delogu & Colafrancesco 1994; Del Popolo & 
Gambera 1996). Each of these effects either increase or de- 
crease the collapse times for the mass shells, of which the 
combined effect is unclear, and such uncertainties should be 
largely absorbed into the ~ 1 — 2 crossing time range as- 
sumed for the relaxation timescale (as is the uncertainty in 
our specific definition of tcross itself). 

We now describe a quantity related to AM/ M which 
does not require following the detailed virialization process 
of the cluster. The mass shell that has just collapsed at the 
epoch of interest, t, reached maximum expansion approxi- 
mately at a time, t/2. For a given peak height, we can invert 
the equation for the time at maximum expansion to obtain 
the initial radius of this shell. The mass enclosed by this 
shell is. 



(21) 



where ri is the initial (proper) radius. Because at maximum 
expansion the shell has not yet crossed the orbits of other 
shells, the enclosed mass M(< Vm) = M{< ri), where r,n 
is the radius of maximum expansion. Let us denote this en- 
closed mass as Mta{t/2) = M{< r^). Analogously, we can 
estimate the mass that has collapsed at the time t — tcross 
by computing Mta{[t — tcross] /2). Hence, an estimate of the 
fractional collapsed mass accreted over the previous crossing 
time is. 



AMta ^ Mta{t/2) - Mta{[t - tcross]/2) 
Mta ^ \ [Mta{t/2) + Mta{[t " tcross]/2)\ ' 



(22) 



which is computed directly from the turn-around times for 
the mass shells and thus does not involve treating the de- 
tailed virialization process. (The crossing time in this case is 
taken to be a constant which is a good approximation for the 
cluster s pro duced by the adiabatic treatment of virialization 



5.2.1 



.) Comparing AMta/ Mta to AM/ M allows us 
to assess the impact of incorporating the virialization pro- 
cess into our calculations. 



4 COSMOLOGICAL MODEL PARAMETERS 

We consider power spectra appropriate for hierarchical clus- 
tering, in particular CDM and power-law P{k). We take the 
linear CDM power spectrum according to BBKS, 



P{k)c 



oc k 



log(l + 2.34g) 
2.Mq 



[1 + 3.89g + {IG.lqf + (5.46g)^ + (e.Tlg)"*] 



-1/2 



(23) 



where we have assumed the scale-invariant form {P{k) oc 
k) for the primordial spectrum. The spectrum is expressed 
in terms of the parameter q = k/(hT), wh ere the "shape 
parameter" is defined to be (Sugiyama 1991:), 



r = f2o/i exp(— S7_B — fisfio 



(24) 



where fis is the density parameter of baryonic mass. We 
adopt a s our standard the value Q.sh^ ~ 0.016 (Copi et 



al. 1995) consistent with Big Bang Nucleosynthesis. How- 



ever, we also examine models with the larger value inferred 
f rom the X-ray gas con tent in clusters, Q, Bh~^ ^ /Q,o = 0.05 
( [White fc Fabian 1995| ). 

We use the shape parameter to set the value of h for the 
CDM models {Ho = lOO/i km s'^ Mpc'^). Viana & Liddle 
( [l996| ) determine 95% confidence limits F « 0.2 — 0.3 using 
the gal axy a utocorrelation function data from Peacock & 



Dodds (1994). We slightly expand this range to F = 0.15 — 
0.35 and construct a linear relation with f2o such that F = 
0.15 for VLo = 0.2 and F = 0.35 for Q.o = 1. This function 
T{^lo,h) - at fixed Q.b - allows us to determine a value of 
/i as a function of Slo consistent with current data (though 
see Peacock 1997) with the slightly expanded range chosen 
to keep the values oi h ~ 0.45 — 0.85 within reasonable 
observational limits (e.g. Fukugita, Hogan, & Peebles 1993). 

To eliminate the divergence of P{k)cDM at small wave- 
lengths we follow the standard procedure and mathemat- 
ically smooth the power spectrum with a gaussian filter, 
P{k)s oc P(fc)cDAf exp(— fc^P/2), where I is the smoothing 
length in comoving coordinates. The mass contained within 
the gaussian fifier is Ms = {2n)^^'^ pbH/y^)'"^ at 2 = 0. 
Since we are interested in mass fiuctuations of the size of 
clusters of galaxies, we choose Ms = 10 h~ Mq which 
corresponds to the mass of a large galaxy at the present 
epoch, since structures with M < M^ individually have 
a negligible contribution to the dynamical evolution of a 
cluster. The smoothing length corresponding to this Ms is 
I = 0.86r2o ^^^/i"^ Mpc. For evaluating AM/M at higher 
redshifts, we approximately account for the smaller masses 
of galaxies by reducing the smoothing mass according to. 
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Ms{z) — Ms(0)/(1 + z), which is just t he self-similar ac- 
cretion law M oc t^^^ ( Bertschinger 1985). (The results we 



obtain for AM/ M in the next section are not overly sensi- 
tive to /.) 

To normalize the smoothed CDM power spectrum we 
first set the value of as = ctm (8/i~^Mpc), the linearly ex- 
trapolated rms mass fluctuation in spheres of radius 8h~^ 
Mpc, to agree with the abundance of X-ray clusters at the 
present day {z = 0), 



0-8 



0-8 



0.52a 



{-0.46+O.lOno) 


(-0.52+0.13no) 



Ao =0 
f2o + Ao 



(25) 



where we have used the best-fit results of Eke, Cole, & Frenk 
(1996). We obtain as for higher reds hifts using the same scal- 
ing relation for the density in §3.1. The number of clusters 



as a function of mass measures the amplitude of fiuctuations 
in the mass on cluster scales (see, e.g., Frenk et al. 1990), 
whereeis we require the amplitude of fluctuations in clusters 
(i.e. integrated over mass). In the peaks formalism, clusters 
are biased tracers of the mass distribution which may be 
accounted for by increasing the mass fluctuation amplitude 
by a bias factor (BBKS; Bardeen et al. 1987), 



'-'clus — 

ao 



(26) 



where {i>} is an average peak height of the model under 
consideration above some threshold ut- The cluster power 
spectrum is then, P{k)cius = b'^i.^^P{k)mass, where P{k)mass 
is the smoothed CDM power spectrum normalized using the 
above relation for as- 

We evaluate bdus in a m anner similar to that described 
in Croft & Efstathiou (|l994|) (also see Efstathiou et al. 1992). 



First, we smooth the mass power spectrum with a gaussian 
of length / — 5V2Q,Q^^^h~^ Mpc and normalize to the above 
erg. This smoothing length gives approp riate cluster masses 



over a w ide range of cluster abundances ( Croft fc Efstathiou 
1994|)7}Ve focus our attention on clusters with masses ex- 
ceeding 3.5 X 10^''/i~^Mq corresponding to an abundance 
of ~ 1 X 10~^ Mpc"^ for typical CDM models (see Figure 
1 of White, Efstathiou, & Frenk 1993). By requiring that 
the number density of peaks in our model match this clus- 
ter abundance, we specify the peak threshold, Vt, and thus 
{{>}. For CDM models with Qq — 0.2 — 1 we obtain bias 
parameters foe! us ~ 1 — 3.5. 

We also consider models with power-law power spectra, 
P{k)pi oc fc", with n ranging from to -2. Because P(k)pi 
also diverges for small wavelengths we smooth it in the same 
manner as P{k)cDM ■ We normalize the smoothed power 
spectrum in the same manner as done for the CDM spec- 
trum and compute appropriate bias factors. This is feasible 
because the as-^o relationship given above is very insensi- 
tive to the shape of the power spectrum (e.g. White et al. 
1993). 

After choosing a power spectrum and specifying the 
remaining model parameters, we construct clusters having 
masses spanning the range (0.35—3) x IQ^^ Mq. However, 
in the peaks formalism the size of an over-density region is 
specified by the peak height, v, rather than the mass. When 
examining clusters at 2 > 0, we use u determined at z = 
e.g. for a 1 X 10^^^~^Mq cluster defined at z = 0, we study 
its progenitor at z > which has the same v but smaller 
mass. Note that whenever we quote cluster masses in the 



next section we refer to the actual mass computed from the 
full spherical model calculation within r = 1.5/i~^ Mpc. 



5 RESULTS 

5.1 Mass profiles at z = 

Although not the focus of our present study, considerable lit- 
erature exists on the theoretical density profiles of clusters. 
Hence, to place our model clusters in context, we briefiy 
summarize the properties of the cluster mass profiles ob- 
tained from the spherical accretion model &t z = Q for a 
CDM power spectrum. A more detailed disc ussion of th e 
density profiles will appear in a separate paper (Buote 1997). 

In Figure |l| we display the density and mass profiles 
of clusters with mass ~ 7 X lO^*/i"^M0 formed in mod- 
els with fio = 0.2, 1. Only the virialized matter is shown 
for the density profiles since the infalling shells contribute 
(formally) infinite density spikes at their locations. We ex- 
press the radii in terms of r2oo, the radius where the mean 
cluster density is 200 times the background value: r2oo = 
1.63 X lO"^(M//i-^M0)^/^ Mpc (e.g. Navarro, Frenk, & 
White 1997). For the clusters shown, r2oo = 1.5/i"^ Mpc. 

Without the infalling material p{r) is virtually identical 
within r-200 for both Q.o — 0.2, 1. For r ^ 0.2r-2oo the virial- 
ized portions of the clust ers approximately follow the self- 



similar profile p ~ r~ ' (Bertschinger 1985) and steepen to 
p ~ r~^'^ at larger radii r ~ (1 — 1.5)/i^ Mpc. When the in- 
falling mass is included the behavior at small radii is largely 
unaffected, but at large radii the flo = 1 profile becomes fiat- 
ter than f^o = 0.2: p ~ r"'^ '* for fio = 1 and p ~ r~^'^ for 
Slo = 0.2 over radn r ~ (0.5 - l)r2oo or r ~ (0.75 - 1.5)/i"^ 
Mpc. For a given Qo the profile shapes of clusters of different 
masses are essentially identical when expressed in terms of 
r200. 

The density and mass profiles appear to agree well with 
those produced in the three-dimensional N-body simulations 
of a standard CDM cosmology. The p ~ j--^/'^ dependence 
in the inner regions giving way to a steeper density profile 
matches the behavior found by Anninos & Norman (199C) 
in their detailed study of the effect of mass resolution on 
cluster density profiles. Moreover, the relative similarity of 
the density profiles for different masses and different values 
of fio agrees with the conclusions of Navarro et al. (1997), 
aside from the core regions {r ^ 0.2r2oo) where Navarro et 
al. find p ~ r~^. 



5.2 AM/M 

5.2.1 CDM: z = (Ao = 0) 

We begin our discussion by examining the results of AM/ M 
at z; = for CDM models with Aq = 0. At present we 
restrict our discussion to the case where the relaxation 
timescale equals tarosa. In Figure |2| we show the radial de- 
pendence of AM/ M for a couple of representative masses 
with rJo = 0.2, 1. The profiles of AM/ M are nearly con- 
stant, though there appears to be a slight (~ 10%) increase 
with radius for low-mass clusters when fJo = 1; i.e. the 
dynamical state, or degree of virialization, is not a strong 
function of radius within 1.5/i~^ Mpc. The slight increase of 
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Figure 1. Density (left) and mass (right) profiles of « 7 X lO^'^h^^ Mq clusters in CDM models with = 0.2 (dashed) and = 1 
(solid). Only the contribution from virialized matter is shown for the density profiles while the mass profiles contain both the virialized 
and inf ailing mass. 
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Figure 2. The radial profile of AM/ Af at z = for CDM 
models with Qq = 0.2, 1. Clusters with mass fa 3.5 X lO^'^h^^MQ 
are shown as filled circles and those of mass ~ 1.4 X lO^^h-^MQ 
by crosses. (These masses are approximately the values within a 
radius 1.5h~^ Mpc.) 



AM/ M with radius is reasonable since the amount of ac- 
creting mass becomes fractionally more important at larger 
radii, and mo dels w ith larger fio have proportionally more 



accretiii g maoo 



* * * The similar radial dependence of AM/ M for different 
masses also indicates that there is no conceptual advantage to 
focusing on a scaled radius like ^200- 



To assess the reliability of AM/ M with respect to 
the simplifications of our model, we compare to the re- 
sults of the N-body simulations of BX. However, BX an- 
alyzed the projection of and focused on the mean val- 
ues of logj^Q Pm/Po {m — 2,3) of their N-body cluster sam- 
ples, where P^/Po = {(4'j^*)')/((*o"')') (see g). Although 
these properties prohibit a rigorous quantitative comparison 
of AM/ M to the BX results, we expect qualitative similar- 
ities. Henceforth we shall compare our results for AM/ M 
to the Pm/Po results obtained by BX with these caveats un- 
derstood. For our comparison at hand, the nearly constant 
radial behavior of AM/ M is similar to Pm/Po, though the 
Pm/Po actually show a modest decrease with radius. This 
decrease in Pm/Po is expected because the properties of the 
multipole expansion dictate that higher order moments in 
the potential decay much more rapidly than the monopole 
with increasing distance from a mass concentration. 

In Figure ^ we display AM/ M as a function of cluster 
mass (i.e. defined within 1.5h~^ Mpc). The dependence of 
AM/ M on cluster mass is significant: for a given radius 
the more massive clusters tend to have smaller AM / M. 
This is a reflection of the shorter relaxation (crossing) times 
for the more massive clusters within a given radius. The 
slightly steeper dependence for the low mass clusters within 
the 1.5/i~^ Mpc aperture just reflects the slight increase with 
radius of AM/ M for Q. = 1 models discussed above (see 
Figure ^ . 

Now we focus our attention on how AM/ M depends on 
SIq. Since observations and N-body simulations give mean 
values averaged over a cluster sample, we average AM/ M 
over the mass function appropriate for our CDM models. Al- 
though the peaks formalism (BBKS) does not provi de a well 
defined mass fu nction, the analytic Press- Schechter ( Press & 



Schcchter 1974) mass function is convenient and a suitable 



approximation in many cases (Bond et al. 1991). For our 
averaging procedure we compute AM/ M for a few clusters 
spanning the mass range (0.35 — 3) x 10^^ h"^ Mq and inter- 
polate AM/ M for arbitrary M over that range. This func- 
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Figure 3. AM/ M evaluated for radii r = 0.75/1"^ Mpc (boxes) 
and r = 1.5h~^ Mpc (crosses) as a function of total cluster mass 
within 1.5h-l Mpc for Uq = 0.2, 1 CDM models at z = 0. 



tion is then averaged over the Press-Schechte r mas s function 
using the approximations in Viana & Liddle ( 1996) for CDM 
models to yield the mass-averaged fractional accreted mass 
denoted by (AM/M). 

In Figure § we display (AA//M) within r = Ih'^ Mpc 
£18 a function of flo {Xo = 0). As expected, (AM/M) in- 
creases with increasing Qo- (AM/ M) = 0.46^0 ''^ when av- 
eraged over the full mass range (0.35 - 3) X 
and (AM/M) ?^ OAOQ,'^-^° when averaged over the smaller 
range (0.7 — 3) x W^^ Mq corresponding to more mas- 
sive clusters. The slight variation of the exponent of Slo as 
a function of mass range is the result of the mass averaging 



- we find that AM/ M cx fio applies essentially over the 
full mass range (0.35 - 3) X IQ^^h'^MQ. Given the resuhs 
of Figures ^ and it is not surprising that there is little 
advantage in going to larger radii as the greatest increase 
is seen within r = 1.5/i~^ Mpc for the least massive clus- 
ters for which we obtain the marginally steeper dependence 
(AM/M) ^ 0.470^''^ averaged over the full mass range. 

1/2 

To determine the relationship between this ~ f^Q 
dependence and the turn-around times of the mass shells 
we compare these results to AMta/ Mta- Setting tcross = 
0.14/ffo (see below) we compute AMta/ Mta as a func- 
tion of collapsed mass aX z — Q. In Figure ^ we dis- 
play AMta/ Mta as a function of Q.o for a collapsed mass 
7x 10"/i"^Mq. In agreement with the behavior of AM/M, 
we find AMta/ Mta oc fio'** which is essentially independent 
of mass over the fuU range (0.35 - 3) x lO^'^ft-^Mg. The 
normalization of AMta/ Mta, however, is ~ 40% smaller: 
AMta /Id ta = Q.29,^o'^^ versus AM/M ^ GA2,nl^ for 
7 X IO^^/i'^Mq. The close agreement between AM/ M and 

AMta/ Mta indicates that the ~ ^^o'^^ dependence is derived 
from the turn-around times of the mass shells. (This depen- 
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Figure 4. The mass-averaged fractional accreted mass evaluated 
for r = lh~^ Mpc at z = for CDM models as a function of 
(Aq = 0). The crosses indicate a mass average over the full range 
(0.35 — 3) X lQ^^h~^ Mq, and the boxes indicate a lower limit of 
7 X 101-*/i-1Mq. Relaxat ion timescales of 1 and 2 crossing times 
are shown. 



deuce is only weakly affected by the relaxation timescale - 
see below.) 

We now compare {AM/ M) to the N-body results of 
BX. In order to make the comparison more reasonable, let 
us consider logio{(AAf/ M)^) within the 0.75/i"^ Mpc aper- 
ture, for which our models give {(AM/M)^) ^ 0.19^0®^ 
over the full mass range. Let us consider the logarithmic 
difference of this quantity between a model with fio = 1 
and ilo < 1 and denote this shift by A. We can then write, 
A — — O.Sllogj^Q f2o- Setting f2o = 0.35 appropriate for BX 
(and Cm ~ 1, see §[2.3|), we obtain A = 0.4, in pretty good 
agreement with A ~ 0.5 obta ined by BX considering that 
they analyzed in projection ]* * *[ This qualitative agree- 
ment with BX is reassuring on two accounts: (1) {AA4/ M), 
which represents the fractional increase in the monopole 
of the potential, clearly does correlate strongly with the 
next few low order moments; (2) our simple spherical accre- 
tion model is able to reasonably reproduce the differences 
in the means of cluster distributions generated by three- 
dimensional N-body simulations. 

We may also compare {AM/ M) to the results of Rich- 
stone et al. (1992) who apply the spherical collapse model 



to clusters of mass lO^^/i~^M0 arising from homogeneous 
density perturbations, which are assumed to form at the 
time equal to twice the turn-around time of the perturba- 
tion. Richstone et al. compute the quantity 5F, the fraction 
of present-day clusters which formed within the last time 



* * * The projection of the dark matter in clusters appears to 
yield smaller variations in Pm/Po than does the projection of 



(see §4 of Tsai & Buote). 
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Figure 5. The fractional "collapsed" mass accreted over the 
previous relaxation time as a function of Co for CDM models 
at 2 = 0. The collapsed mass is approximately 7 X W^'^h~^ Mq 
enclosed within a radius 1.5/i~^ Mpc. The solid lines indicate the 
Ao = models, and dashed lines indicate f2o + Aq = f models. 
Relaxation timescales of 1 and 2 crossing times are shown where 

tcrosa = O.IA/Ho. 



interval St, as a function of f2o and Aq. Consulting Figure 3 
of Richstone et al. we see that 5F ~ Q.o for substructure sur- 



vival titne dt — U.L/Uo (,a timescale similar to our itc 



(A similiii' dSJi^itld^itiCS 6{ 5F on f2o is sdkii in Fig'urS 13 6f 
Lacey & Cole 1993.) 

Although 5F is computed from the collapse times of 
density perturbations, it depends on these collapse times 
through the error function; i.e. {AM/ M) and 5F are quali- 
tatively different in how they depend on the collapse times of 
shells and thus f2o- In fact, we can repeat the exercise above 
for computing the mean shift between models with f2o = 1 
and f^o < 1 by replacing (AM/ M) with 5F in which case 
we obtain A = —2 logjQ fio- For Q.q = 0.35, A = 0.9 which is 
almost twice the mean shift found by BX. The significantly 
better agreement of the A predicted by (AM/ M) over 5F 
with the results of BX demonstrates that a much closer rela- 
tionship exists between {AM/ M) and Pm/Po than between 
SF and Pm/Po as would be expected. 

The crossing times for the (Ao — 0) CDM models are not 
a strong function of mass or Qq'- tcross = (0.12 — 0.17)/-ffo 
with a typical value tcroas — O.IA/Ho = 1.4 x 10^h~^ yr. 
(The largest value corresponds to a mass 3.5 x IO^^/i-^Mq 
for f^o = 1.) We show in Figure ^ the results for (AM/ M) 
for the case where the relaxation timescale is set to 2t cross ■ 
Although the values of (AM/ M) have approximately dou- 
bled, the power-law dependence on f2o remains virtually un- 
changed: (AM/M) ^ l.Qml^^ when averaged over the full 
mass range (0.35 - 3.5) x lO^'^/i'^M©, and (AM/M) ^ 
0.94ng'*^ when averaged over the smaller range (0.7 — 
3) X lQ^^h~'^MQ. This behavior is mirrored by AM^/Mta 



which maintains its fio"'® dependence for 2taross (Figure ^. 
(Note this same type of fio-dependence also applies for a 
relaxation timescale of 0.5tcross.) 

A noticeable weakening of the dependence on fio oc- 
curs when the relaxation timescale is increased to iit cross- 
(AAf/M) ?^ l.iml^^ when averaged over the full mass 
range (0.35-3.5) x IO^'^/i'^Mq, and (AM/M) ^ 1.6m[]''^ 
when averaged over the smaller range (0.7-3)x10^''/i-^Mq. 
For a cluster that accretes an increasingly larger fraction 
of its current mass during the relaxation time, AM/ M ap- 
proaches a value of 2 independent of f2o- Because the smaller 
masses tend to have larger AM/ M, their dependence on Q.q 
weakens more rapidly with redshift than the more massive 
clusters. 

AMta/ Mta is slower to respond to this effect. For 
3tcross, AMta/Mta = 1.2ir2o'"^ for a mass 3.5 x 



lO^'^h'^MQ. This becomes l.SOfio ''^ at Atcross and approx- 
imately 2.09^0^'^ at 5tcross. The departure of AMta/T^^ta 
from AM/ M indicates that the dynamical evolution de- 
pends more on the virialized structure of the cluster for large 
relaxation times, or equivalently for large r redshifts. (This 
effect is discussed in more detail in 55.2.3.) 



For our current purposes the trend of (AM / M) w ith f2o 



is of principal importance, not the actual value . Recall 

that in this paper we do not intend to compare (AM/ M) 
directly to observations (§ 2.3), but rather we expect the next 
few moment ratios of the gravitational potential, which can 
be observed, to be strongly correlated with (AM/ M) (see 
f^. Since 1-2 crossing times should effectively cover the rea- 
sonable range of relaxation timescales for potential fluctua- 
tions characterized by the lower order moments, the utility of 
(AM/ M) for studying the dependence of cluster morphol- 
ogy is hardly affected by the uncertainty in the relaxation 



timescale. This behavior contrasts wit h that of (Rich, 
stone et al. 199^ ) and related measures (Lacey & Cole 199J ) 



whose value is to be directly compared to the "frequency of 
substructure" . 

Changing the value of (equation |2^) has virtually 
no effect on_AM/M. For crs = 0.52,1 (Hq = 1) we flnd 
that AM/ M varies by 1%. This similarity occurs be- 
cause the amplitude of P(k)cius arising from a higher erg 
is mostly compensated for by a reduced cluster bias factor, 
bcius- This insensitivity to cts = 0.5 — 1 is also found in the 
mean logj^Q Pm/Po distributions of BX. 

Finally, we investigated the effe ct of using Qb appropr i- 
ate for the gas fractions in clusters (White & Fabian 1995). 



In our models this only has the effect of slightly changing 
Ho determined for a particular value of flo (see definition of 
F in As a result, (AM/ M) is essentially unaffected - 
we find results in excellent agreement with those obtained 
above with the BBN Qb- Of course, it would be better to 
test the effects of a greater baryon fraction by incorporat- 
ing the gravitational dynamics of a dissipational component 
into the model. 
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Figure 6. AM/ M versus radius for models with Qq = 0.2, 
Ao = 0.8 (boxes) and Co = 1 (crosses). 
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Figure 7. (AM/ M) versus Hq for f2o + Ao = 1 models (boxes) 
and Ao = models (crosses). The relaxation timescale is Itcross- 



5.2.2 CDM: fio + Ao = 1 at z = 

We now turn our attention to the zero-curvature models, 
f^o + Ao = 1. In a low-density universe with zero-curvature 
only recently has the efltect of a cosmological constant sig- 



* * * This trend is nu ite in,s ensitive to the slope of the power 
spectrum — see end of jj5.2.4| 



nificantly imprinted itself on the cosmic dynamics; i.e. the 
transition from f2o = 1 to f^o < 1 occurs considerably later 
and much more rapidly than in the Ao = case. Hence, 
although the shells which collapsed at the present epoch 
turned around at times similar to those in Ao = mod- 
els, the previously collapsed shells turned around with times 
governed by a universe with fio ^ 1- 

As a result, the density profiles of the zero-curvature 
models have essentially the same slope as the f2o = 1 mod- 
els over radii ~ (0.5 — l)r2oo (see §5.1). However, the pro- 



file of the zero-curvature models is flatter at smaller radii, 
r ^ 0.2r2oo, indicating that the mean density within the cen- 
tral regions of the clusters is considerably smaller in zero- 
curvature models, consi stent with the arguments given in 
§4 of Richstone et al. (1992). Because of this decrease in 
density, the half-mass radii are larger, the velocity disper- 
sions smaller, and thus the crossing times (equation |l^) are 
considerably larger in the low-density zero-curvature models 
than in their Ao = counterparts: tcross ~ (0.20 ~ 0.23) / Hq 
over the full mass range for the Qq = 0.2 and Ao — 0.8 model, 
a factor of ~ 1.7 larger than the open Qq ~ 0.2 model. In 
fact, with respect to f^o = 1, tcross is ~ 40% larger for the 
zero-curvature ilo = 0.2 model but ~ 15% smaller for the 
open Jlo = 0.2 model. This considerably stronger and quali- 
tatively different dependence of AM/ M on crossing time as 
a function of ilo causes the low-density zero-curvature mod- 
els to have qualitatively different behavior than their Ao = 
counterparts. 

We display in Figure ^ the radial proflle of AM/ M for 
the Slo ~ 0.2, Ao — 0.8 model and the ^Iq — 1 model for a 
couple of representative masses. For M ^ 7 x 10 h~ Mq, 
the profile of the zero-curvature model is essentially flat but 
almost of the same magnitude as the f2o = 1 model. As 
expected, the differences between the models becomes most 
noticeable at larger radii because the accretion of the most 
recent shells is suppressed in the zero-curvature model; i.e. 
the recently collapsed shells approximately behave as in the 
Qo ~ 0.2, Aq = case. 

The Qo = 0.2 zero-curvature radial profile of AM/ M 
slightly decreases with radius for M — 3.5 x 10^ h~ Mq and 
actually exceeds the Qo — 1 values for r < Ih"^ Mpc. This 
decreasing profile indicates that shells with turn-around 
radii ^ 2h~^ Mpc dominate the contribution to AM/ M 
over the previous crossing time - since the crossing time 
is significantly larger than for 57o = 1, the zero-curvature 
model accretes shells at higher redshift which turned-around 
at smaller radii. 

In Figure]^ we display {AM/ M) versus Q,o for the zero- 
curvature models. In order to accentuate the dependence on 
Qo we focus on the largest aperture considered, r = 1.5/i~^ 
Mpc. We find that these models behave very nearly as 
Qo = 1 and depend only weakly on Qo: {AM/ M) ~ Qo'^ . 
This behavior is qualitatively consistent with the linear 
growth of density fluctuations: the ratio of linear growth 
factors for the Qo = 0.2, Aq = 0.8 model to f2o = 1 is 0.70 
compared to a ratio of 0.33 for the open Qo — 0.2 model to 
Qo — 1. That is, from linear growth we would expect the 
amount of accreted mass in clusters in the Qq = 0.2 zero- 
curvature model to be more similar to the Qq = 1 model 
than to the open Qq = 0.2 model (e.g. Peebles 1984). How- 
ever, the N-body simulations of BX indicate that the means 
of the logjQ Pm/Po distributions of their zero-curvature mod- 
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els, though intermediate with the Slo = 0.35, 1 cases they ex- 
amined, nevertheless leaned more closely to the f2o = 0.35 
open model for radii (0.4 - 1.2)/i~^ Mpc. 

We emphasize that this weak f2o-dependence arises 
from the increased crossing times for the low-density zero- 
curvature models. In Figure ^ we display AAfta / Mta for the 
zero- curvature models with tcross = 0.14/_H'o appropriate for 
the Ao = models. In this case AAfta /Mta oc Oq'^^ which is 
only slightly weaker than the f2o-dependence for the Ao = 
models. This ~ fio '*^ dependence applies to {AM/ M) if we 
fix tcross ~ 0.14/_ffo for the zero-curvature models. (This 
behavior is also found for SF by Richstone et al. who also 
used a constant relaxation timescale for all models.) 

Our procedure of setting the relaxation timescale equal 
to the crossing time is only an approximation. Because the 
crossing times in the zero-curvature models vary by almost 
a factor of 3 more than in the open models over Qq ~ 0.2—1, 
the accuracy associated with this procedure is considerably 
more important for studying the f2o-dependence of AM /M 
in the zero-curvature models. 

In addition, other factors contribute to a computational 
uncertainty. Because the increase in crossing time is due pri- 
marily to the flatter density profile at small radii, the error 
due to the number of shells becomes more important. Also, 
issues like the smoothing length of the power spectrum be- 
come more important at radii ^ 0.1r2oo. Hence, the detailed 
features in the profile of {AM/ M) vs Q,q in Figure ^ (e.g. 
the r^o ~ 0.7 values exceeding f2o = 1) are probably due to 
the uncertainty in the relaxation timescale. Further study of 
the relaxation timescale for the zero-curvature models using 
three dimensional N-body simulations would certainly be of 
value here. 

Nevertheless, it is clear that the low-density zero- 
curvature models do not mimic the corresponding open mod- 
els and thus cluster morpholo gies d o indeed distinguish be- 
tween these two cases at 2 = Ol* * 1. We find that the differ- 
ence between the models arises from the detailed virialized 
structure of the clusters (i.e. crossing times), which is con- 
sistent with the expectations of Lahav et al. (1991) that a 
cosmological constant, "could only have changed the density 
profile by its efi'ect on the non-linear motion of shells." 



5.2.3 CDM: z>Q 

Let us now consider the evolution of AM/ M with red- 
shift in the context of the CDM model. For compactness we 
restrict our discussion to a relaxation timescale of Itcross 
and fis equal to the BBN value. (The effect of different 
relaxation timescales is mentioned at the end of this sec- 
tion.) Since the cluster population becomes less massive 
with increasing z, we analyze AM/ M averaged over masses 
[(0.7-3)7(1-1-2)] X IQ^^h-^Msun, where the [l + z) factor is 
intended to account for the reduction in mass according to 



*** BX found that the means of the log^gPrn/Po distribu- 
tions of 39 simulated clusters for their Oq = 0.35 open and 
zero-curvature models, though systematically larger for the zero- 
curvature models as indicated by the non-overlapping error bars 
derived from bootstrap re-sampling, were formally consistent in 
terms of the Student's-t test. It appears that a larger cluster sam- 
ple is required to formally differentiate the mean s of t he two mod- 
els. This has been confirmed by Thomas et al. (1997). 



the self-similar accretion law ( Bertschinger 198E). We focus 



on a slightly higher mass range because a feasible observa- 
tional sample either from lensing or X-rays will be biased 
towards more massive clusters. The qualitative results pre- 
sented in this section remain if the lower limit on the mass 
average is halved. 

In Figure ^ we display the radial profile of {AM/ M) 
at 2 = 0.4, 0.8 for models with f^o = 0.2, 1 (Aq = 0) and 
fio = 0.2, Ao = 0.8. Over redshifts 2 ~ - 0.4 the behav- 
ior of {AM/ M) with radius does not change appreciably. 
The shape of the open fio = 0.2 profile remains constant 
over this interval, but the f2o = 1 profile shifts from being 
slightly increasing with radius (Figure ^) to slightly decreas- 
ing. As discussed in § 5.2.2 for the zero-curvature models, the 



decreasing profile indicates that as 2 increases the clusters 
in the f2o = 1 models are growing faster in their inner re- 
gions; i.e. the accreted mass is dominated by shells with 
turn-around radii ^ 2h~^ Mpc. This contrasts with the sit- 
uation at 2 = for the Ao = models where the accretion 
of infalling shells contributes almost entirely to large radii, 
which is consistent with the behavior of the exact self-similar 



solution (Bertschinger 1985). 

The shape of the fractional accreted mass profile of the 
zero-curvature fio = 0.2 model is similar to that of the cor- 
responding open model, but the normalization is similar to 
Slo = 1 as we found for the 2 = case (with the s ame c aveats 
regarding the crossing time efi'ects discussed in § 5.2.2| ). The 
profile of the zero-curvature Q.o = 0.2 model does not de- 
crease like the fio = 1 model because the shells which af- 
fect most the outer radii are still determined by f2o ~ 0.2. 
That is, for the zero-curvature models the redshift indicat- 
ing when the term involving the cosmological constant in 
the Friedmann equation dominates the matter term is given 



by, (1 -)- Ztrans) 



fin 



1, after which linear growth is 



suppressed (e.g. §13 Peebles 1980). For 57o = 0.2, we have 



Ztr 



0.6. 



Indeed, proceeding to higher redshift, 2 = 0.8, we see 
in Figure ^ that the profile of the zero-curvature f^o = 0.2 
model now decreases with radius indicating that the mat- 
ter term is beginning to dominate the Friedmann equation. 
The {AM/ M) radial profile for f2o = 1 steepens while that 
for the open Q,o — 0.2 model remains essentially flat. For 
redshifts ^ 1 the profiles of the zero-curvature Qo = 0.2 
model and fio = 1 model have very similar shapes. Even the 
open model starts to steepen, though from the argument 
from linear theory we would not expect the transition to oc- 
cur until approximately ztrans = 3. We emphasize, however, 
that the precise redshifts of these occurrences depend on the 
relaxation timescale. For longer relaxation timescales, these 
properties occur at lower redshifts and vice versa. 

For the higher redshifts it is clear that the smaller aper- 
ture sizes optimize differences in Qq (Ao = 0). In Figure ^ 
we show {AM/ M) computed within a radius 0.5/i~^ Mpc 
as a function of fio for 2 = 0.4, 0.8. Although the value of 
{AM/ M) increases with increasing redshift indicating that 
the cluster morphologies are more disturbed, its sensitiv- 
ity to Q.0 decreases significantly with redshift: {AM / M) « 
O.eOSTig-^^ for 2 = 0.4 and {AM/M) ^ 0.87n|^-^^ for 2 = 0.8. 
(For comparison, for r = 1.5h~^ Mpc we obtain {AM/ M) ~ 
O.em'^-^^ for 2 = 0.4 and {AM/M) ^ 0.75n[i °^ for 2 = 0.8.) 
The weakening trend with redshift is even observed as early 
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Figure 8. The radial profile of (AA// M) as a function of redshift for CDM models Qq 
(boxes), and Qq = 1 (crosses). See text for details of the mass average. 



: 0.2, Ao = (filled circles), Ho = 0.2, Aq = 0.8 



in the lh~^ Mpc 



as z = 0.2 where {AM/ M) 
aperture. 

This decrease in the sensitivity to fio arises because of 
two effects. First, as z increases the dynamics of the uni- 
verse approaches that of the Einstein-de Sitter case regard- 
less of the present value of ilo and Aq. Second, the relaxation 
timescale becomes an increasingly larger fraction of the age 
of the universe, and thus {AM/ M) —> 2 regardless of Qq 
and Aq. 

However, the decline in sensitivity to fio is more rapid 
for AM/ M than for AMta/T^ta- We find that the relation- 

1 /2 

ship AMta/ Mta ~ f2g holds from 2 = until z « 2 at 
which point AMta/ Mta ~ 2 and the Slo-dependence van- 
ishes. Since AMta/ Mta depends only on the turn-around 
times of shells and not the detailed virialized cluster struc- 
ture, its scales more similarly to what is expected from linear 
theory: as remarked above, for the open flo — 0.2 model the 
transition to the Einstein-de Sitter phase occurs approxi- 
mately at ztrans = 3. However, this redshift is not reached 
because the crossing time is essentially the age of the uni- 
verse at z « 2. Clearly, taking into account the detailed 
virialized structure of the clusters is important for higher 
redshifts. 



5.2.4 Power law 

Now we switch our focus to models having power-law power 
spectra. In Figure we display AM/ M vs spectral index 
n at z = for clusters of mass 7x lO"/i"^M0; other masses 
have very similar profiles. For the moment we concentrate 
on models with 57o = 1 and have a relaxation time of Itcross. 
The value of AM/M for n « -1.3 agrees with the = 1 
CDM value as expected for clusters. The fractional accreted 
mass increases as n decreases and varies more over the range 
n = 0, —2 (r^o = 1) than it does as a function of fio for CDM 
models over the range fio = 0.2, 1; AAf/M « 0.13(-n)^''^ -f 
0.15. 

This n-dependence may be understood in terms of the 
peaks formalism as follows. Peaks of a given height in n = 
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Figure 9. {AM/M) evaluated for r = 0.5/i"^ Mpc at z = 
0.4,0.8 for CDM models as a function of f2o (Aq = 0). See text 
for details of the mass averaging. 



models are much more isolated than peaks of similar height 
in n = —2 models; i.e. a cluster of mass ~ is an 

exceedingly rare high peak ~ 50 in n = models whereas 
such a cluster is a "normal" ~ 3 peak for n = —2 models. 
Models with smaller n have more power on small scales and 
hence more possible mass to accrete at later times, which is 
to be regulated by the value of JIq. Thus, regardless of the 
value of fio = 0.2 — 1, the huge peak in the n = model has 
a smaller value of AM/ M than the normal peak in n = —2 
model. 

However, real clusters are not i/ ~ 50 peaks. Here we 
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Figure 10. AAf/ M for cluster of mass 7 X IO^-^/i-^Mq eval- 
uated for r = lh~^ Mpc at 2 = for power-law models with 
Q.Q = 1 as a function of spectral index n. 



Figure 11. AM/ M vs Qq for power-law models with n = — 1, —2 
for clusters of mass 7 X lO^'^h-'^MQ evaluated for r = l/t~i Mpc 
at 2 = 0. 



experience the principal flaw in the peaks formalism in that 
merging is not taken into account (although see Bond & 
Myers 1996), which particularly affects models with n > — 1 
because in this case the collapse dynamics of a peak is not 
dominated by its initial density distribution ( Bernardcau 
1994|)~]We have argued (see §p|) that the effects of merg- 
ing should not affect the mean of AM/ M when averaged 
over a cluster sample which is most relevant for studying 
Qo- Rather, the fluctuations induced in AM/ M by mergers 
necessarily affects the variance. 

To understand this qualitatively let us envision clus- 
ters formed in n = 0, —2 models considering the effects of 
mergers. In the n — model, since there is little mass on 
smaller scales for a cluster of reasonable peak height to ac- 
crete, AM/ M will tend to be small for many clusters. How- 
ever, when a cluster merges with another of comparable size 
(it will be similar in size since there is little small-scale mass) 
AM/ M will be comparatively large. In comparison, clusters 
in the n = — 2 model have significant small-scale mass to ac- 
crete so AM/ M will not vary as much from cluster to clus- 
ter as in the n = model; i.e. when averaged over a sample 
of n = clusters, AM / M will have a broader distribution 
than in the n = —2 model. 

The three dimensional N-body simulations of BX 
support this qualitative picture: the distributions of 
log^QPm/fb for the n = Q model indeed have significantly 
larger variances than for n — —2. (But the means of the 
distributions are largely the same.) Hence, for studying the 
effect of the slope of the power spectrum on the relationship 
between AM/ M and Q,o merging must be properly taken 
into account. Since our spherical model does not account 
for the effects of mergers, it is imperative to keep the shape 
of the power spectrum approximately fixed when comparing 
models with different fio as we have done in the previous sec- 



tions. It is also imperative to insure that the effective n on 
cluster scales is less than —1 so that the collapse dynamics 
i s determined most ly by the initial peak density distribution 
(Bcrnardeau 1994). 

Although it is important when comparing AM/ M for 
models with different Slo to insure that n is approximately 
fixed, the value of n is not important. In Figure hll we show 
AM/ M as a function of fio for power-law models with 
n = —1,-2. Although the values of the n = —2 models 
are approximately double those of the n = —1 models, the 
dependence on Qq is similar in each: AM/ M ~ O.Slf^o'^'^ 
for n = ~l and AM/M^ O.gSO.'^ *^ for n = -2. 



6 DISCUSSION AND CONCLUSIONS 

We have constructed a simple, intuitive model to study the 
dependence of substructure in clusters on Qq, Ao, and z. 
We characterize the importance of substructure and non- 
ellipsoidal features in the mass of a cluster, or equivalently 
the dynamical state, in terms of, AM/ M, the fractional 
amount of mass accreted over the previous relaxation time 
within a radius r, where AM is the mass increase and M 
is the average mass within r over the relaxation time. This 
fractional accreted mass is the monopole term in the multi- 
pole expansion of the fractional increase in the gravitational 
potential due to matter interior to r. Since, by hypothe- 
sis, AM/ M determines the importance of substructure and 
non-ellipsoidal features, it must strongly correlate with the 
next few multipole terms which are approximately the ra- 
tios, $j"7$™«. 

These multipole ratios, or rather their projected coun- 
terparts, P,n/Po = {(*m')^)/{(*o"*)^), are directly observ- 
able and share the same type of correlation with AM/ M 
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(see §2.2). With N-body simulations the correlations be- 



tween AM/ M and Pm/Po can be explicitly quantified and 
thus allow a measurement of flo to be obtained directly from 
cluster observations with gravitational lensing and X-rays. 
In this paper we have examined the differences of AM/ A/ 
in models with different Qq (and P{k)) which are related to 
the differences in Pm/Po without requiring knowledge of the 
actual correlation. (The reasonable agreement found for the 
mean shift between 57o = 1 and f2o = 0.35 for the Pm/Po 
of the N-body clusters in BX and AM/ M pr ovides quan- 
titative justification for this procedure - see §5.2.1.) That 



is, we predict how the low order Pm/Po should behave with 
Qq, Xo, and z. In this way we can acquire intuition of how 
clusters evolve dynamically in different cosmologies, and, as 
we discuss below, we can identify the best ways to constrain 
no, Ao with much greater ease than via computationally ex- 
pensive three dimensional N-body simulations. We compute 
AM/ M using the spherical accreti on mo del in the adiabatic 
approximation of Ryden & Gunn (L987). 

Our model extends and improves upon previous stud- 
ies in several respects. First, we compute complete one- 
dimensional mass distributions of clusters using the spheri- 
cal accretion model, where the initial density profiles of the 
clusters are given by the peaks formalism (BBKS). The re- 
laxation time, taken to be a multiple of the crossing time, 
is computed individually for each cluster. In contrast, pre- 
vious studies (Richstone et al. 1992; Kauffmann & White 
1993; Lacey & Cole 1993) did not compute mass distribu- 
tions of individual clusters and assumed a constant relax- 
ation timescale for all clusters. Moreover, Richstone et al. 
assumed clusters arise from homogeneous density fluctua- 
tions on one particular mass scale. 

The principal difference, however, between our work 
and the previous studies mentioned above lies in how they 
relate to observations. The previous studies predict the frac- 
tion of clusters that are expected to show evidence of sub- 
structure. In terms of observations this is generally referred 
to as the "frequency of substructure" and interpreted as the 
fraction of observed clusters that exhibit multiple density 
peaks. Because these previous theoretical studies do not in- 
dicate the type of substructure or morphological features to 
be seen in clusters, it is unclear how to compare their pre- 
dictions with observations; i.e. what sizes of multiple density 
peaks are required before a cluster is to be considered un- 
formed? This ambiguity is resolved in our model because 
AM/ M is expected to correlate strongly with for 
small I and with ^'™*/*o"' for small m. These multipole 
ratios quantify the cluster morphology explicitly. 

Because the previous studies do not specify the type 
of substructure, the relaxation timescale (or rather the sur- 
vival time of substructure) is very uncertain - with some 
s tudies indicating uncer tainty up to ~ 1 — 10 crossing times 
(Nakamura et al. 1995). Since the predicted "frequency of 
substructure" is highly sensitive to the substructure survival 
time, the predictions made by the previous studies are too 
uncertain to be useful (Kauffmann & White 1993; Lacey & 
Cole 1993) . In our model, the relevant relaxation timescale is 
realistically confined to ~ 1 — 2 crossing times since AM/ M 
is concerned only with the low order moments in the cluster 
potential (i.e. large potential fluctuations containing ^ 20% 
of the total cluster mass). Moreover, since it is only the ex- 
pected correlation of AM/ M with the low order moments 



that we require, and not the actual value, the uncertainty 
associated with ~ 1 — 2 crossin g times is unimportant for 
our present study of Oq (§5.2.1). 

The results of our study suggest the following proce- 
dure to optimize constraints on flo and Ao from analysis 
of cluster morphologies quantified in terms of Pm/Po- The 
present epoch exhibits the strongest dependence of cluster 



morphologies on Qq with {AM/ M) ~ n, 



1/2 



the higher 



mass clusters M ^ 7 x 10^* ft ^-^0 ^.re slightly preferred for 
this purpose. At z = Q, {AM/ M) depends only weakly on 
radius. However, Pm/Po declines with radius because the 
higher order moments give way to the monopole term as 
the radius is increased. Hence, any r ^ Ih^^ would seem 
best in this light. (Although not too small to invali date the 
connection between the 2-D and 3-D monopoles - §2.2.) 

Although the dependence of {AM / M) weakens with 
increasing redshift, analyzing high-redshift data separately 
and combining with z = Q data significantly improves the 
ability to constrain models. First, simply obtaining inde- 
pendent cluster samples at medium {z ~ 0.4) and high 
[z ~ 0.8) redshifts and combining with z = data in- 

1/2 

creases the dependence on f2o from ~ at 2: = 



to 



X 



X f2o ~ 0,0. Second, the evolution 
of {AM/ M) with z depends on Q,o as well: the ratio of 
{AM/ M) for z = 0.8 and z = Q models goes approxi- 
mately as SIq indicating that the morphologies of clus- 
ters in a low-density universe undergo more rapid evolution 
than if flo = 1 thereby providing an additional constraint. 
Finally, the radial dependence of {AM/ M) is very different 
for f^o = 0.2, 1 models over radii (0.5— 1.5)/i~^ Mpc, and we 
would expect that Pm/Po would decrease with radius con- 
siderably more rapidly for z ^ 0.5 if Slo = 1 than if Oo = 0.2 
with or without a cosmological constant. 

Because the dependence of AM/ M on r^o weakens with 
increasing redshift. X-ray observations of clusters have a 
clear advantage over gravitational lensing observations since 
the latter are essentially restricted to 2 <; 0.15. Moreover, 
the Pm/Po computed on the X-ray surface brightness distri- 
butions of clusters have a larger dynamic range than those 
computed on the projected mass distributions (§4 of Tsai & 
Buote 1996), which probably trans late s to a greater sensi- 
tivity for probing fio (see end of §2^). For z <; 0.2, both 
X-ray and lensing observations are well suited for probing 
the evolution of cluster morphologies with redshift. How- 
ever, the theoretical comparison for lensing data only re- 
quires dissipationless simulations which is a definite advan- 
tage over the more computationally expensive dissipational 
simulations required for the comparison to X-ray data. 

Although we believe at present that further refine- 
ments to our model are not justified given the uncertainty 
in the relaxation time and the qualitative nature of the 
relationship between AM/ M and the low order multi- 
pole ratios, there are areas which could use some atten- 
tion. A study quantifying the relationship between AM/ M, 
(($-')2)i/2/(($-*)2)i/2^ and Pm/Po usiug N-body Simula- 
tions would help clarify predictions of this model (i.e. de- 
termine Cm and dm - see §2.2), and will enable a direct 



comparison to observations of cluster surface density maps 
obtained from gravitational lensing and X-ray data. If such 
a study can also significantly limit the range of relaxation 
times for low order potential fluctuations, then investigation 
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of the secondary effects such as non-sphericity and angular 
momentum may be warranted. 

To illustrate the potential for such studies, let us com- 
pare A M/ M to the Pm/Po computed in §4 of Tsai & 
Buote ( 1996| ) for the dark matter (i.e. p as opposed to 
p^) of clusters formed in a f2o = 1 CDM simulation. The 
even Pm/Po are nearly proportional to each other such that 
P2/P0 « IOOP4/-P0, and the mean value of P2/P0 is ~ 10"''. 
If we assume the same scaling applies to the monopole, i.e. 
{AM/Uf ^ IOOP2/P0, then we estimate AM/M « 0.10. 
This value is in qual itativ e agreement with the value of « 0.4 
we determined in 



5.2.1 



Since the magnitude of AM/M 
scales approximately linearly with relaxation time, this sug- 
gests that the appropriate relaxation timescale for the low 
order multipoles is < Itcross = 0.14/_ffo. 

Hence, cluster morphologies would appear to be an at- 
tractive means to probe Slo and Ao. Theoretically, cluster 
morphologies are related to Qq and Ao in a conceptually 
straightforward way, are sensit ive to Ao at z = 0, and are 
very insensitive to Qb (5:5.2.1) and the cluster bias factor 
unlike some other indicators for Qq (Dekel et al. 1996). Ob- 
servationally, cluster morphologies are easy to quantify in 
terms of Pm/Po (Buote & Tsai 1996). The challenge lies in 
acquiring lensing and X-ray data of large samples (<: 50) of 
clusters at various redshifts for comparison to large, high- 
resolution N-body simulations of lensing maps and hydro- 
dynamic simulations of X-ray images. 

Alternatively, if N-body simulations confirm that the 
correlations between AM/ M and the low order Pm/Po are 
tight and mostly independent of Q.o, then the N-body sim- 
ulations will only be needed to initially calibrate this cor- 
relation. In this case subsequent comparison of AM/ M to 
observations can be achieved rapidly using the spherical ac- 
cretion model and allow Q,o and Ao to be determined with 
considerably less computational expenditure than required 
for the standard method. 
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APPENDIX A: 
CASE J2o + Ao = 



SPHERICAL MODEL FOR 
1 



Here we give the formulae for the radius and time at max- 
imum expansion of spherical shells evolving in a universe 
where f2o 4- Aq = 1 (see §[j.2|). Related expressions may be 
found in Richstone et al. (19921) for the case of homogeneous 



density perturbations, an d in E ke et al. (1996), who adopt 



the equations of Peebles ( 1984 ) which assume a particular 
initial condition, a = r at z ^ 00, where a is the expansion 
factor and r is the proper radius of the mass shell. We give 
expressions explicitly defined for an arbitrary initial epoch. 

We begin with the energy integral for a spherical per- 
turbation in a universe with non-zero cosmological constant 
A = 3//2a„ 



s=im^^M(lI)_lA.^ (Al) 

2\dtJ r 6 ' ^ ' 

which is conserved before shell-crossing. The mass enclosed 
within proper radius r is. 



M{<r) = ^p,{U){l + 5,)rl 



(A2) 



where pb{ti) is the density of the background universe and 
the proper radius at the initial time ti which has expanded 
to radius r at a time t < tm- The energy evaluated at ti is. 



(A3) 



where Ki — Hfri/2 is the initial kinetic energy and we 
have neglected peculiar velocities and set Ai = 1 — f^i. At 
maximum expansion the energy is all potential. 



^ {-K,n, [1 + 5,]) + (—) {~K, [1 -n,]), (A4) 



where rm is the proper radius at maximum expansion. Since 
the energy is conserved we may set Em = Ei and obtain an 
equation for Vm/ri^ 

^ — I / T".__ \ " n ■ I T... \ 

= 0, (A5) 



- 1 _ i — \ + 1 

1 + 5i \ril 1 + 5i / 

solutioi 
(3 + 3) 



which has the standard cubic solution, 

1/2 



yfi y^-/:^ - 1 



where, 



cosO = 



3^/^ i+I. / n-i-i y^^ 

2 - 1 I 5,: ) ' 



(A6) 



(A7) 



for < < 7r/2. Finally, by setting the energy before max- 
imum expansion equal to the initial energy we obtain an 
equation for the time as a function of r. The time corre- 
sponding to maximum expansion is thus, 

tm = U^Hl^Q.-^''^(\+~& 



'1/2. 



dRVR 



n- 



1 + S, 



-1/2 



(A8) 
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